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Abstract. A multivariable generalisation of the Bessel polynomials is intro- 
duced and studied. In particular, we deduce their series expansion in Jack poly- 
nomials, a limit transition from multivariable Jacobi polynomials, a sequence 
of algebraically independent eigenoperators, Pieri type recurrence relations, 
and certain orthogonality properties. We also show that these multivariable 
Bessel polynomials provide a (finite) set of eigenfunctions of the hyperbolic 
Sutherland model with external Morse potential. 



1. Introduction 

The Bessel polynomials first appeared in a paper by Bochner [Boc29] as the 
polynomial solutions of the differential equation 

2 d2 y , I , u\ d y (a i\ 
x + (ax + bj — = n(n + a — l)y, 

ax z ax 

where a and b are (real) parameters, and n a non-negative integer. The first 
systematic and detailed study of their properties was later undertaken by Krall 
and Frink [KF49]. For further details and references see e.g. the book by Gross- 
wald [Gro78]. In the present paper we introduce and study a multivariable gen- 
eralisation of the Bessel polynomials as eigenfunctions of the partial differential 
operator 

(i) ^=£^+£(^+^+ 2 ^ 



dx 2 ^— ' dxi t-r 1 Xi — Xj dx 



i—l 1 i—1 i^j 



with k a (real) parameter. In particular, we obtain their series expansion in Jack 
polynomials, a limit transition from Jacobi polynomials associated with the root 
system BC n , n algebraically independent eigenoperators, Pieri type recurrence re- 
lations, and, as further discussed below, certain orthogonality properties. We note 
that by setting n = 1 in {I} we indeed recover the eigenoperator of the one- variable 
Bessel polynomials. As in the one- variable case, the dependence on the parameter 
b can be removed by substituting (6xi, . . . , bx n ) for (x-y, . . . , x n ). Following the 
standard convention we will accordingly set 6=2. 

The differential operator |T]) is closely related to the Schrodinger operator of the 
so-called hyperbolic Sutherland model with external Morse potential. The results 
we obtain on the multivariable Bessel polynomials are thus relevant also for this 
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quantum many-body model. In order to make this relation precise we let 

n 

* (xi, . . . ,x n ) = H xf- 2V2 e-^ Y[( Xl - Xj ) K . 



i—l i<j 

As a direct computation shows, conjugation of D B by this function an d a 
change of coordinates to (21, . . . , z n ) = (e Xl , . . . , e Xn ), yields the Schrodinger type 
operator 

( 2) -E^+E («--(»-)«-) 



with 



2k(k - 1) V ^ 

4-*. 4sinh 2 Uzi - Zi) 



#o = -^(n 2 - 1) - -i \ ^ " n{n - 1) - V - j n; 



3 v , 2 4 

see e.g. Hallnas and Langman [HL07] for a proof of this fact. The Schrodinger 
type operator ^ was studied by Inozemtsev and Meshcheryakov [IM86]. They 
determined the discrete part of its spectrum, and showed that the corresponding 
eigenfunctions are given by symmetric polynomials. We also mention that the 
(hyperbolic) Sutherland model, obtained from ((2|) by eliminating the external Morse 
potential, was introduced and studied by Sutherland [Sut71] . 

It is clear from (J5J that any eigenfunction of D B yields an eigenfunction of H 
upon multiplication by the function ^o- The physically relevant eigenfunctions are 
those contained in the Hilbert space L 2 (M. n , dzi, . . . dz n ). From this point of view 
the eigenfunctions of D B contained in L 2 (R™ , | "P^ol 2 , ■ ■ ■ dx n ) are of particular 

interest. In Section 7 we will show that all multivariable Bessel polynomials up to 
a given degree are contained in this latter Hilbert space if and only if the param- 
eters a and k satisfy a simple inequality. We then prove that these multivariable 
Bessel polynomials form an orthogonal system with respect to the Hilbert space 
inner product. In addition, we deduce an explicit expression for the corresponding 
squared norms. In effect, this provides a (finite) set of normalised eigenfunctions of 
H, representing bound states in the corresponding quantum many-body system. 

Multivariable generalisations of the classical orthogonal (Hermite, Laguerre and 
Jacobi) polynomials, of the same type as the multivariable Bessel polynomials con- 
sidered here, have been extensively studied in the literature. We mention, in par- 
ticular, the closely related work of Lassalle [Las91a, Las91b, Las91c] and Macdon- 
ald [Mac], as well as Baker and Forrester [BF97] and van Diejen [vD97,vD99]. For 
further references see e.g. the book by Dunkl and Xu [DX01]. As is well known in the 
one-variable case, the Bessel polynomials share a number of properties with these 
classical orthogonal polynomials, e.g., the fact that they solve a second order ordi- 
nary differential equation. The present study of multivariable Bessel polynomials 
thus provide a natural complement to the extensive literature alluded to above. We 
also mention that multivariable generalisations of the Bessel polynomials, different 
from the one considered here, have been introduced and studied by Exton [Ext93], 
as well as by Pathan and Bin-Saad [PBOO]. 
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We conclude this introduction by a brief overview of the paper. In Section 2 
we briefly review some basic facts on Jack- and multivariable Jacobi polynomi- 
als, as well as generalised hypergeometric series, which will be used in the paper. 
In Section 3 we give a precise definition of our multivariable Bessel polynomials 
and deduce their expansion in Jack polynomials. In Section 4 we establish a limit 
transition from multivariable Jacobi polynomials associated with the root system 
BC n . This limit transition is then used in Section 5 to obtain n algebraically in- 
dependent eigenoperators of these multivariable Bessel polynomials. In Section 6 
we continue to apply this limit transition to obtain Pieri type recurrence relations. 
As already mentioned above, the orthogonality of the multivariable Bessel polyno- 
mials is discussed in Section 7. In Section 8 we present some further remarks on 
the relation between the multivariable Bessel polynomials and the eigenfunctions 
of the hyperbolic Sutherland model with external Morse potential. In this section 
we also discuss the problem of constructing a (moment) functional with respect to 
which all multivariable Bessel polynomials would be orthogonal. 

2. Preliminaries 

In this section we recall some basic facts on Jack polynomials and multivariable 
Jacobi polynomials associated with the root system BC n , as well as on generalised 
hypergeometric series, which we will make use of in later parts of the paper. 

2.1. Symmetric functions and Jack polynomials. We recall that a partition 
A = (Ai, A2, . . .) is a finite sequence of non-negative integers such that 

Ai > A 2 > • • • . 

The number of non-zero parts is commonly referred to as the length of A, and 
denoted by ^(A). For simplicity of exposition we will not distinguish two partitions 
that differ only by a string of zeros at the end. The weight of a partition A is the 
sum 

|A| = Ai + A 2 H 

of its parts. The diagram of A is the set of points € 1? such that 1 < j < A^. 
By reflection in the main diagonal the conjugate partition A' is obtained. We let 
be the unit vector in Z™ defined by {e{)j = Sij, and write 

A w = A + a, A w = A — ej 

for i = 1, . . . , n. On the set of partitions the dominance order can be defined by 

fi < A \fi\ = |A| and fj,\ + ■ ■ ■ + fn < Ai + • • • + A,, Vi = 1, . . . , n. 

We stress that two partitions A and /x are comparable only if they have the same 
weight, i.e., if \/j,\ = |A|. The monomial symmetric polynomials m\ are given by 

m x {%i i • ■ • i x n ) = ^ ' ■ ■ ■ x n n , 

q£S„A 

where S n X denotes the orbit of A under the natural action of the permutation group 
S n of n objects. These monomials span the space 

A„ C C[xi, ...,x n ] 
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of symmetric polynomials in the variables x = (xi, . . . , x n ) with complex cocffi- 
cients. For any non-negative integer m we let 

m 
fe=0 

where Ajj C A„ is the subspace of homogeneous polynomials of degree k. For further 
details on partitions and symmetric functions see e.g. Chapter I in Macdonald 
[Mac95]. 

In order to proceed we require the differential operators 

d_ 

dxi 



i=i 
and 

d 2 „ v — ^ x? 



OX: . , . Xi Xj Qx% 

1=1 1 J 

for I = 0, 1 and k = 1,2, respectively. We recall that the (monic) Jack polynomials 
P\ can be defined as the unique symmetric polynomials that satisfy the following 
two conditions: 

(1) P x =m x + X^<a u Xfl m^ for some coefficients u Xfl , 

(2) P\ is an eigenfunction of E\ and D2; 

see e.g. Macdonald [Mac95] or Stanley [Sta89]. To any two partitions A and ^ a 
generalised binomial coefficient ( ) can be associated via the expansion 

Pa(xi + l,...,x n + l) -r-^ f\\ P\(xi, . . . ,x„) 



= E 



We recall that the fact that the sum extends only over partitions \i C A is non-trivial. 
To the best of my knowledge, this fact was first established by Kaneko [Kan93]. 
A short and elegant proof, exploiting a relation with the so-called shifted Jack 
polynomials, was later given by Okounkov and Olshanski [0097]. 

For future reference we state the known action of the operators Eg and Dk , as 
well as the power sum p\(x) = x\ + ■ ■ ■ + x n , on the Jack polynomials: 

f3al F PX{X) - Vf A ^ PAw(x) 

(3b) E\P\ = \\\P X , 

/q V n P\(X} ( A \(\ 1 J. ( 'W P \^) 

3c Di—— — r= > , (Ai — 1 + K(n — ij) ' . , 

1 ' Px{i n ) y W)/ ^A (l) (l") 

(3d) D 2 P x = d x P x , rf A = ^A J (A J -l + 2 K (n-z)), 

(3e) 

where we have made use of the notation 

h x = Y\Y\(\ l -j + K{\' J -i) + i), 
i=i j=i 
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and where the sums are over all i — 1, . . . ,n such that in (|3a|) and (|3c[> A(,-) is a 
partition and in (|5e|) \^ is a partition. We remark that (|3b[) and {3d} are readily 
inferred from the definition of the Jack polynomials as common eigenfunctions of E\ 
and Z?2 , and the action of these operators on the monomials m\ . Equation (|3a[) can 
be deduced from the definition of the binomial coefficients (see e.g. Kaneko [Kan93]), 
and (pc|) can then be obtained using the fact that 2D\ = [Eq, D%\. Equation ([3ej) . on 
the other hand, is more difficult to prove. As far as I know, it was first established 
by Kaneko [Kan93]; see also Okounkov and Olshanski [0097]. 



2.2. Generalised hypergeometric series. Kaneko [Kan93], Koranyi [Kor91] and 
Macdonald [Mac] introduced and studied multivariable hypergeometric series as- 
sociated with Jack polynomials. In order to recall their definition we require the 
following generalisation of the Pochhammer symbol: 



1{X) 



Y[[a-K(i-l)} Xi , 



where the ordinary Pochhammer symbol [a] m is given by [a]o = 1 and [a] m — 
a(a + 1) • • • (a + m — 1) for m > 0. The generalised hypergeometric series can now 
be defined by its series expansion in Jack polynomials: 



(4) p F q (ot 1 ,...,a p ;Pu---,0q;x]K) - ^ 



i4 kL '[4 k) ^(^) 



A LM1JA Ih'qlX 

For p = 2 and q = 1 it is known that the resulting generalised hypergeometric series 
2F\{a.\, ct2] Pi', x; k) satisfies the diferential equation 

(5) (£>i - D 2 )F + - k{h - l))E F - (a x + a 2 + 1 - n{n - l))EiF = n ai a 2 F; 
see e.g. Kaneko [Kan93] or Yan [Yan92]. 



2.3. Jacobi polynomials associated with the root system BC n . Jacobi poly- 
nomials associated with root systems have been extensively studied in the literature. 
We mention, in particular, the work of Vretare [Vre84], the first to introduce such 
polynomials, and that of Debiard [Deb88], Heckman and Opdam [H087], as well as 
that of Beerends and Opdam [B093]. We will make use of the Jacobi polynomials 
associated with the root systems BC n , henceforth referred to as the BC n Jacobi 
polynomials. These polynomials are even and permutation invariant trigonometric 
polynomials, and can be defined as eigenfunctions of the differential operator 

dBC = E §j + E ( fci coth \* + 2k * coth z >) ^ 

+ h y: (coth \{ Zi - Zj) (± - A) + coth \{ Zl + Zj) (J- + 
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with ki, ki and k% (real) parameters. As discussed by Beerends and Opdam [B093], 
by a change of variables to t j = — sinh Zi/2 this operator can be written as 

D BC = - !)^2 (h + k 2 + - ~ + 2k 2 + l)tA 

i=l i i=l 

k(U-l) d 



2k 3 J2 



U — tj dt 



i¥=3 

It is thus clear from formulae ([3|) that 

D BC P x {t{z)) = e BC P x (t(z)) + £ c x ,P,(t(z)) 

for some coefficients c> M , and with 

e BC =d x + (k 1 +2k 2 + l)\\\. 

This leads us to the definition of the BC n Jacobi polynomials P BC as the unique 
polynomials satisfying the following two conditions: 

(1) P^ c {z) = {-±) W P\{t{z)) + Y,^ cX uxnPx{t{z)) for some coefficients u x „, 

(2) P^ c is an eigenfunction of D BC . 

We remark that it is common to define the BC n Jacobi polynomials in terms of the 
monomials m x (e Zl + e~ Zl , . . . , e z ™ + e~ Zn ); see e.g. Beerends and Opdam [B093]. 
However, using the triangular structure in the expansion of the Jack polynomials 
in the symmetric monomials m x it is easy to see that this leads to an equivalent 
definition. As will become apparent below, for our purposes the present definition 
is somewhat more convenient. We also mention that the leading coefficient of the 
BC n Jacobi polynomials P x c has been set to (— 4)' A ' in order to ensure that they 
are monic when expanded in the monomials specified above. 



3. Definition and series expansion in Jack polynomials 

In this section we will first give a precise definition of our multivariable Bessel 
polynomials, and then proceed to deduce and study their expansion in Jack poly- 
nomials. We start by observing that ([3|) implies that the operator D B acts on the 
Jack polynomials as follows: 

with 

(7) e x = d x + a\X\ 

(where we recall that 6 = 2). This implies, in particular, that D B maps each Jack 
polynomial P x to a linear combination of Jack polynomials P u labelled by partitions 
/i C A. It follows that D B have eigenfunctions of the form 

Y x (x; k) = ^ u^P^x; k) 
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for some coefficients u\^. indeed, using ((6]) it is readily verified that this is the case 
if and only if the coefficients satisify the recurrence relation 

(8) {e x ~ e„)u x „ = 2 ^ p ^ ^ j « v<) 

for all partitions /i C A. If e M 7^ e> for all such partitions, then the coefficients 
with C A are uniquely determined by ((8|) once the value of u\\ is fixed. This 
leads us to our definition of multivariable Bessel polynomials. 

Definition 3.1. Fix the values of the parameters a and n, and let A = (Ai, . . . , A„) 
be a partition such that e M ^ e\ for all partitions /i C A. We then define Y\ as the 
(unique) symmetric polynomial such that 

(1) Y\ = P x + J2^c\ U *» P » f or some coefficients ux^, 

(2) Y\ is an eigenfunction of D B . 

It is clear from the discussion above that the level of non-degeneracy required of 
the eigenvalues e\ is essential in order for the multivariable Bessel polynomials Y\ 
to be well-defined. As we now show, this condition is generically satisfied. 

Proposition 3.1. Let A = (Ai, . . . , X n ) be a partition. Then, for generic values of 
the parameters a and K, ^ e\ for all partitions /1 C A. In particular, this is the 
case if k > and a satisfies the condition 

(9) o<-2(|A| + /c(n-l)) + l. 

Proof. We fix a partition /i C A. Since e\ — e M is a polynomial in k and a, and not 
identically zero, we have that e M ^ e\ for (k, a) in an open dense subset of R 2 (or 
even C 2 ). We observe that A is of the form 

n 

A = /i + 2_. n i e n n>i € N. 
i=l 

It follows that 

n 

e\ - e M = 2J ni (2{jii + n(n — i)) — l + a). 

i=l 

Under the stated conditions on k and a the right hand side is clearly non-zero. □ 

Throughout the remainder of the paper we will assume that the parameters 
a and k are such that the Bessel polynomials in question are well-defined. We 
mention, already at this point, that in our discussion of orthogonality in Section 
7 the relevant parameter values will be precisely those for which k > and a 
satisfies the condition ^j. Proposition 13.11 will then ensure that the multivariable 
Bessel polynomials under consideration indeed are well-defined. With that said we 
proceed to solve the recurrence relation ([8]), and thereby to obtain an expression 
for the coefficients in the series expansion of the multivariable Bessel polynomial in 
Jack polynomials. By iterating ([5]) we find that the coefficients ux^ are given by 

where the sum is over all chains of partitions 

T: A = (°'AD (1) Ad ...D «A = /i 
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such that r — |A| — and |^ l_1 'A| — |^'A| = 1 for all i = l,...,r, i.e., over 
all standard tableux T of shape A//i; see e.g. Section 1.1 in Macdonald [Mac95]. 
Inserting the expression ([7]) for the eigenvalues of the Bessel polynomials we thus 
arrive at the following: 

Proposition 3.2. The multivariable Bessel polynomials have a series expansion 



with coefficients 

(10) U\f2 

T 

where the latter sum is over all standard tableux T of shape A//i. 

Remark 3.1. For the multivariable classical orthogonal polynomials this type of 
series expansions are well known. To the best of my knowledge, such a series ex- 
pansion was first obtained by Constantino [Con66] for the multivariable Laguerre 
polynomials corresponding to k = 1/2. We recall that for this particular param- 
eter value the Jack polynomials specialise to the so-called Zonal polynomials; see 
e.g. Macdonald [Mac95]. James and Constantine [JC74] later obtained the analo- 
gous series expansion for the multivariable Jacobi polynomials with k = 1/2 or 1. 
Series expansions for general values of k can be found in Baker and Forrester [BF97] , 
Beerends and Opdam [B093], Lassalle [Las91a,Las91b,Las91c], as well as in Mac- 
donald [Mac]. 

The generalised binomial coefficients are in general rather complicated ob- 
jects; see e.g. the remark in Section 3 of Okounkov and Olshanski [0097]. However, 
the particular binomial coefficients appearing in Proposition 13.21 can be written 
down explicitly using the fact that 

,„, er)=(^' + ^»-'»n \-_\ + :: + t;:o l> - 



This latter formula can be deduced from a comparison of (|3e[) with the correspond- 
ing Pieri formula in Macdonald's book [Mac95] (Equation 6.24 (iv) in Chapter VI); 
see Lassalle [Las98] for further details. We recall that Stanley [Sta89] (see Theorem 
5.4) obtained the specialisation of the integral form 

<?(A) A, 

Ja = ^h x P x , h x = Ul[(X^j + K (X' 3 ~ i + 1)), 
i=i j=i 

of the Jack polynomials at x = (l n ). Combining this specialisation formula with 
the appropriate limit of Proposition 3 in Lassalle [Las98] (multiply by 1 — t, set 
q = t a and let t — > 1) it is readily verified that 



P\W (1") _ „ \ + k(ti - i + 1) -q Aj - Xj + 1 + n(j 



(12) 



P A (l n ) A, + k(£(X) + L} X l -X J + l + K (j 
i _ A, - Xj + n(j + 



n 



Xi-Xj + Kij-i) 
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Using these two formulae we can thus write down an explicit expression for the 
coefficients u\^. This enables us to study the dependence of the multivariable Bessel 
polynomials on the parameters o and k. In particular, since the coefficients in the 
expansion of the integral form of the Jack polynomials in terms of the monomials 
m\ are polynomials in l/n (see e.g. Knop and Sahi [KS97] and references therein), 
we can prove the following: 

Corollary 3.1. The multivariable Bessel polynomials are rational functions of a 
and k, and are thus for generic values analytic functions of these parameters. In 
particular, this is the case if k > and a satisfies the inequality ^ . 

Proof. It is immediate from the discussion above that the coefficients are ra- 
tional functions of a and n. The first part of the statement can thus be established 
via an argument similar to that used to prove Proposition 13.11 The second part 
of the statement follows from Proposition [3Tj inspection of (fTTj) and (112|) . and the 
fact that the Jack polynomials P\ are analytic in k for k > 0. □ 

In general, it seems hard to deduce a closed expression for the coefficients u\^. 
However, if A = (k n ) for some non- negative integer k the situation simplifies con- 
siderably, as the multivariable Bessel polynomials then are given by a terminating 
generalised hypergeometric series. More precisely, we have the following: 

Proposition 3.3. For any non-negative integer k, 

F(fen)(a;; a, n) = Ck(a, K,) 2 F (~k, k + a - 1 - n(n - 1); — at/2; k), 

where 

n , , fyfc")(«) 

y ' r 1.1 M \h _i_ „ _ i _ _ im(«) 



-k]^ n) [k + a-l-K{n-l)]\ kn) 

ittii 

k]{ K) {a 2 }^ P,(x) 



Proof. Since (—k)\ = unless A C (fc n ), we obtain upon setting p = 2, q = 1 and 
a% = — k in (0J that 



Ac(fc") IPiJa 

The fact that the Jack polynomials P\ are homogeneous of degree |A| implies that 
2 F (-k,a 2 ;-x/2) = lim 2 Fx(-k,a 2 ; f3 1 ;-(3 1 x/2). 

01— »oo 

Substituting —(3\x/2 for a; in (J5]), and then taking the limit /3j — > oo, we thus find 
that 2 Fq(— k, a 2 ; —x/2) satisfies the differential equation 

D 2 F + (a 2 + 1 - k - K(n - l))EiF + 2E Q F = nka 2 F. 

It follows that 2 Fo(—k, k + a — 1 + n{n — 1); —x/2) is an eigefunction of D B . We 
now obtain the statement by comparing leading coefficients of this generalised hy- 
pergeometric series and Yr \. □ 

4. A LIMIT TRANSITION FROM BC n JACOBI POLYNOMIALS 

In this section we establish a limit transition from the BC n Jacobi polynomials 
to the multivariable Bessel polynomials. We will make use of the representation 

d bc _ e BC 

(13) P? C [z) = II e B C _J c (-4) |A| iM*(*)) 
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for the BC n Jacobi polynomials P BC . It is easily verified that the polynomials 
defined by the right hand side of (fl3|) indeed satisfy the two defining properties 
of the BC n Jacobi polynomials. First of all, it is immediate from ([3]) that the 
right hand side of (|I3p is triangular in the Jack polynomials P\ . In addition, since 
the operator W^xiD® — e BC ) annihilates the subspace spanned by the Jack 
polynomials P n with /.i C A, it is also clear that the right hand side of (fT5|) is an 
eigenfunction of D BC . We mention that this type of representations have been used 
by Stokman and Koornwinder [SK97], as well as van Diejen [vD99], to study limit 
transitions between different families of multivariable orthogonal polynomials. 

It is clear that also the multivariable Bessel polynomials Y\ have a representation 
of the form (fT3|) . obtained by substituting D B for D BC , e B for e BC , and P\{x) for 
(— 4)! A 'P\(i(z)). As we will show below, there exist a limit transition from the BC n 
Jacobi polynomials P BC which, in effect, performs these substitutions. Using the 
fact that 

e 2z + l e z I 

coth z = —= = 

e 2z — 1 e z — 1 e z + 1 

we can rewrite the operator D BC as follows: 



=i 1 i=i 



* 3 £ 

i<j 



h - 1 e z * + 1 J J dzi 
q q \ e zi+zj + 1 / 9 9 



dzi dzj J e z i+ z i — 1 V dzi dzj 



If we set 

(14) h = ~ {a - 1 + 2e e ) , k 2 = j (a - 1 - 2e e ) , k 3 = k, 

and substitute (z\ + e, . . . , z n + e) for (21, . . . , z n ), then we obtain the limit 

lim D BC = V ^ + V (a - 1 + 2e~ z *) + 2 K V 

i— 1 * z— 1 

which is precisely the eigenoperator D B of the multivariable Bessel polynomials 
written in the coordinates (zi, . . . , z n ) = (lnxi, . . . , lnx„). Moreover, we have 

ef C = ef , lim e~ c sinh 2 ^-t_f = £_. 

By applying the limit in question to (fT5|) we thus obtain the following limit transi- 
tion from the BC n Jacobi polynomials to the multivariable Bessel polynomials: 

Proposition 4.1. With the parameters (ki, k 2 , k 3 ) given by (|14[) we have that 

lim e-^P BC ( Zl +e,...,z n + e)= Y x (e z \ . . . , e z "). 

e — >oo 

5. Higher order eigenoperators 

In this section we prove that the multivariable Bessel polynomials are common 
eigenfunctions of n algebraically independent partial differential operators. We will 
obtain this result by applying the limit transition in Proposition 14.11 to a sequence 
of eigenoperators for the BC n Jacobi polynomials constructed by Heckman [Hec9I]. 
We remark that Heckman obtained such eigenoperators for the Jacobi polynomials 
associated with any (integral) root system. In addition, we study the structure of 
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the commutative algebra generated by these higher order eigenoperators of the mul- 
tivariable Bessel polynomials, and establish an Harish-Chandra type isomorphism 
onto the algebra of even and permutation invariant polynomials in n variables. 

We start by briefly reviewing Heckman's construction, specialised to the root 
system BC n . In doing so we will essentially follow the exposition in Section 2.3 in 
Matsuo [Mat92]. The first step consists of inductively defining sequences of partial 
differential operators Df^, wnere d £ N and i — 1, . . . , n, by setting Dq ± { — 1 
and 

(15) D B d c ±l = ± AdBO i± . ± (h C oth \ Zl + k 2 cothz^ (D B ^ ±t - £>|5 liTi ) 
for <i > 0. If we sum over ±i then we obtain the partial differential operators 



n BC _ I V V n BC 



£ = ±1 j=l 

where we have inserted a factor 1 /2 to simplify certain formulae below. We mention 
that the corresponding sum of the operators D B d+i ±% ' s zer0 f° r all ^ G N. This is 
a direct consequence of the fact that Df^ — (— \) d D B which is easily inferred 
from the definition above. It was shown by Heckman [Hec91] that the BC n Jacobi 
polynomials are common eigenfunctions of the operators D BC , 

We proceed to apply the limit transition to the multivariable Bessel polynomials. 
Accordingly, we let the parameters {k\,k2,k^) be given by (|14[) and substitute 
(zi +e, . . . , z n + e) for (z\, . . . , z n ) in (|15p . In analogy with the discussion in Section 
4, we then obtain the limit 

(16) D B ±l = ± Xl JL D B _ h±z ± ~(a - 1 + 2xr x ) (D B _ h±l - D B _ 1>Ti ) 

± ? E ^ W-i* - ^-i.*) ± f E W-i.* - ^-r,«) 

with -D,f ± .j = 1 and, as before, (x\, . . . , £„) = (e Zl , . . . , e z "). The corresponding- 
limit of the operators D B is simply given by 

1 " 

Dfl = T E E ®2d,ei- 

e=±l i=l 

It is a straightforward exercise to verify that Df = D B , i.e., that the second order 
eigenoperator D B of the multivariable Bessel polynomials is the first element in 
this sequence of operators. Referring to Proposition 14. 1 1 it is now easy to prove the 
following: 

Theorem 5.1. The multivariable Bessel polynomials Y\ satisfy a system of differ- 
ential equations 

D B Y x = e B (X)Y x , d=l,2,..., 
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with eigenvalues of the form 

e B d {\) = \\ d + --- + \ 2 n d + ld. 

(where l.d. stands for terms of lower degree). Moreover, the differential operators 
D B with d = 1, . . . , n are algebraically independent. 

Proof. The fact that the multivariable Bessel polynomials Y\ are eigenfunctions 
of the operators D B is a direct consequence of Proposition 14.11 and the discussion 
above. We observe that, in the order of partial differential operators Xjd/dxj with 
j = 1, . . . , n, the leading term of D B is 




This clearly implies the stated form of the eigenvalues e B . Since the power sums 
Pd{x) = xf+- ■ -+xf l with d = 1, . . . , n are algebraically independent (see e.g. Section 
1.2 in Macdonald [Mac95]), we thus obtain the statement. □ 

Corollary 5.1. The differential operators D B with d = 1, ...,n preserve A n , and 
are thus well-defined as linear operators on A n . As such, they commute with each 
other. 

Proof. We fix a non-negative integer m, and consider the linear subspace A^™\ We 
assume that k > and that a satisfies the condition © . According to Proposition 
13.11 the multivariable Bessel polynomials Y\ with |A| < m are then well-defined. 
Moreover, their triangular structure in the Jack polynomials implies that they span 
A^ m ; c.f. Section 1.6 in Macdonald [Mac95]. Under the assumed conditions on k and 
a, the differential operators D B are thus well-defined, and pairwise commutative, 
as linear operators on A^ m . Since the coefficients of each D B are analytic in k and 
a, this fact immediately extends to arbitrary (complex) parameter values. Since m 
was chosen arbitrarily, the statement thus follows. □ 

Theorem 15.11 provides rather little information on the eigenvalues e B . In partic- 
ular, it reveals nothing of their dependence on the parameters a and k. Indeed, in 
general, it seems that this dependence is rather complicated. However, by apply- 
ing the limit transition in Proposition 14. II to results obtained by Heckman [Hec91] 
(see Theorem 3.11) and Heckman and Opdam [H087] (see Proposition 2.9) on the 
Jacobi polynomials associated with root systems, it is relatively easy to establish 
the existence of eigenoperators of the multivariable Bessel polynomials which have 
eigenvalues of a very simple form. On the other hand, these eigenoperators do not 
seem to have a simple explicit expression. 

In order to make these remarks precise we introduce the commutative algebra 
of differential operators 9 n = C[L>f , . . . , D%]. We let p B = (pf , . . . , p B ) be the 
vector given by 

(17) pf = n(n + (a -l)/2, i = l,...,n. 

To each D £ Sl n we associate a polynomial "/(D) £ C[a?i, . . . , x n ] by requiring that 

(18) DY x = 1 {D){\ + p B )Y x 

for all partitions A = (Ai, . . . , A n ). Since j(D) (A + p B ) is given by a polynomial in 
A, it is clear that 7(D) is indeed uniquely defined by (fT8|) . We will use the notation 
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C[xi, . . . , Xn) 1 " 2 * Sn for the algebra of even and permutation symmetric polynomials. 
Then we have the following: 



Lemma 5.1. For all D 6 



7(D) e C[x, ■■ :;y ' s " 



Proof. We will obtain the statement by applying Proposition 14. II to the analogous 
result for the BC n Jacobi polynomials. We start by recalling this latter result. For 
that we require the vector p BC — (p B c , . . . , p BC ) defined by 

p BC = k 3 (n-i) + (fei + 2fe 2 )/2, i = l,...,n. 

For each differential operator D € C[Df c , . . . , D B c ] we let j BC (D) € C[xi, ...,x n ] 
be defined by 

(19) DP B c = i BC (D)(X + p BC )P? c 

for all partitions A = (Ai, . . . , A„). Proposition 2.9 in Heckman and Opdam [H087] 
then implies that r y BC (D) 6 C[xi, . . . , £ n ] z 2 y,s "-. We observe that with the param- 
eters (^1,^2,^3) given by (fT4"]) we have p BC = p B . We now obtain the statement 
as a direct consequence of Proposition 14.11 and the discussion preceding Theorem 
EU □ 

We can thus define a so-called Harish- Chandra algebra homomorphism 

7 :^„^C[xi,...,x„] z ^ 5 " 

by letting each D £ 2# n be mapped to 7(D). Just as in the case of the BC n Jacobi 
polynomials, the algebra homomorphism 7 thus defined is in fact an isomorphism 
of algebras. That it is injective is immediate from the fact that the operators D B 
are algebraically independent, and that it is surjective is a direct consequence of 
the following analogue of the first part of Theorem 3.11 in Heckman [Hec91]: 

Proposition 5.1. For any polynomial p £ C[xi, . . . , XnY" 2 yiSn there exists a dif- 
ferential operator D p £ 2> n such that 



D p Y x =p(X + p B )Y ) 



A 



for all partitions A = (Ai, . . . , A„). 



Proof. We recall that the power sums p r (xf, . . . ,x%) with r = 1, . . . , n generate 
C[xi,..., x n ) Z2 ™ . By induction in the degree of the polynomials p the statement 
thus follows from Theorem 15. II and Lemma I5TT1 □ 

Remark 5.1. It is interesting to note that the differential operators D p can be 
constructed in terms of the so-called Cherednik operators. To this end we let 
Di denote the Cherednik operator associated with the root system BC n and the 
unit vector e^; see e.g. Cherednik [Che91] or Opdam [Op95] for its definition. We 
furthermore let D BC denote the restriction of p(Di, . . . , D n ) to C[a;i, . . . , £„] z 2 ™ Sn . 
Theorem 2.12 in Opdam [Op95] then implies that D p can be obtained by applying 
the limit transition in Proposition 14. II to D BC . The differential operators D B arise 
in the same way, but with Heckman's 'global' Dunkl operators (see Definition 2.2 
in [Hec91]) substituted for the Cherednik operators. The reason why the operators 
D p lack the simple recursive structure of the operators D BC can now be understood 
as follows: for a Cherednik operator Di there exist w £ ZJJ x S n such that wDi ^ 
D w q\w. In the case of the 'global' Dunkl operators, on the other hand, equality 
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always holds, and it is precisely this property which is responsible for the simple 
recursive structure of the differential operators D±!?; see e.g. Section 2.3 in Matsuo 
[Mat92]. 

As observed above, we thus have the following: 

Corollary 5.2. 7 is an isomorphism of algebras. 

In our discussion of orthogonality in Section 7 we shall need the following: 

Corollary 5.3. Fix a positive integer m, and let A = (Ai,...,A n ) and p = 
(fix, . . . , p n ) be two partitions such that \X\,\p\ < m and X ^ p. Assume that 
K > and that a satisfies the condition 

a < -2(m + n{n - 1)) + 1. 

Then there exist D € & n such that 7(D) (A + p B ) ^ j(D)(p + p B ). 

Proof. By Proposition 15.11 it is sufficient to prove that there exist a polynomial 
p £ C[xx, ■ • • , xj z 2 * s " such that p(A + p B ) ^ p(p + p B ). We observe that 

(A + p B ) n <...<(A + p B ) 1 <0, 

and similarly for p. Since the polynomials in C[xi, . . . , Xn] 1 " 2 y>Srl separate points in 
the open subset of R™ defined by the inequality x n < ■ ■ ■ < xx < 0, the statement 
thus follows. □ 

Remark 5.2. We note that there are values of k and a, and partitions A and p, for 
which there exist no D G Sl n that separates A and p. For example, with k = —2 
and a — 1 this is the case for A = 3(5 and p = 5, where 5 = (n — 1, n — 2, . . . , 0). 
However, it is clear that such examples are rather rare. 

6. Recurrence relations 

It is well known that the one-variable Bessel polynomials y n satisfy the three- 
term recurrence relation 

[J2n + o)(2n + a - 2)| + a - 2j (2n + a - l)y n (x) 

= (n+a — l)(2n + a - 2)y n+ x{x) - n(2n + a)y n -x{x) 
when normalised such that 

Vn{x) = + ^ 1 n Y (n) (x); 

see e.g. Krall and Frink [KF49]. In this section we will obtain a natural multi- 
variable analogue of this recurrence relation: the explicit expansion of all products 
between an elementary symmetric polynomial and a multivariable Bessel polyno- 
mial in terms of the multivariable Bessel polynomials themselves. We mention 
that this type of recurrence relations are often referred to as Pieri type formu- 
lae. Since the elementary symmetric functions generate A n (see e.g. Section 1.2 
in Macdonald [Mac95]), this enables us to obtain the expansion in multivariable 
Bessel polynomials of the product between any symmetric polynomial and a mul- 
tivariable Bessel polynomial. In deducing these recurrence relations we will start 
with the simplest case, corresponding to the elementary symmetric polynomial of 
degree one. We will sketch a proof of this recurrence relation, based on the or- 
thogonality of the multivariable Bessel polynomials, as deduced in Section 7, and 
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their explicit series expansion in Jack polynomials, as obtained in Proposition 13. 21 
We will then proceed to treat the general case using recurrence relations obtained 
by van Diejen [vD99] for the BC n Jacobi polynomials, and the limit transition to 
the multivariable Bessel polynomials established in Proposition 14.11 In order to 
simplify a comparison with his results on the BC n Jacobi case, we will to a large 
extent make use of the same notation as van Diejen. 

As will become apparent below, it will be convenient to employ a normalisation 
of the multivariable Bessel polynomials different from the monic Y\ . This particular 
normalisation has a simple expression in terms of the functions 



i<3 



where 



Af (z) = Y[d* ± ( Zi - Zj )dl ± {z % + « j )II<±(«i). 

i=l 

t(±k + z + s±) 



<±(*) = 



T(z + s±) ' 
r(±(a- l)/2 + z + s±) 

r(2z + s± ) 



with s + = and s_ = 1. Indeed, we can now define 'renormalised' multivariable 
Bessel polynomials by 

, A B (o B ) 

(20) Yx = 2-M +y > Y x , 

A B (p B + X) 

where the vector p B is given by (|17p . In the one- variable case it is straightforward to 
verify that this particular normalisation of the Bessel polynomials coincide with the 
normalisation employed by Krall and Frink [KF49], i.e., Yi n \(x) = y n (x). It is also 
interesting to note that the normalisation in question has a simple characterisation 
in terms of a specialisation formula. 

Proposition 6.1. The multivariable Bessel polynomials Y\ are such that 

(21) n(o n ) = i, 

i.e., the constant term ofY\ equals one. 

In order to prove the statement we will make use of the limit transition from the 
BC n Jacobi polynomials, and a known specialisation formula for these polynomials. 
To recall this formula we let A BC denote the functions obtained from upon 
substituting 

zbci \ _ r(±fc 3 + z + s ± ) 

dv ' ±{z) - T(z + s ± ) ' 

;bc _ r(±(fci + 2fc 2 )/2 + z + s±)r(±(fci + l)/2 + z + a±) 

°^ ±(Z) - r(2z + s± ) 

for d B ± and d B ± , respectively, where, as before, s+ = and s_ = 1. It has been 
shown by Opdam [Op89] (see Corollary 5.2 for R = BC n ) that 

>BC(„BC 



pBC (Q n \ = o2|A| A + (P +^) 

A 1 ; A BC ( P B C) 
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with the vector p BC = (pf c , ■ ■ ■ , p BG ) given by 

pf c = k 3 (n-i) + (fei + 2fc 2 )/2, i = 1, . . . , n. 

We remark that Opdam formulated his result in terms of a function c which (in the 
case of the root system BC n ) is defined somewhat differently from A+. However, 
using the duplication formula T(z)T(z + 1/2) = T(2z)^/ : k /2 2z_1 for the Gamma 
function it is a straightforward exercise to verify that the two formulations are 
equivalent. 

Proof of Provosition \6.1\ We let the values of the parameters (k\, ^3) be speci- 
fied by HH). By definition, 



for some coefficients u\^. Since lim^oo e e t{z + e) = — e z /4, we have 

lim e- £ l A lpf c (z + e) = V(-4)"W ( lim e -*(|A|-|/i|)„ A ^ PJ e z ). 

It follows from Proposition 14. II that 

F A (0") = lim e~ e l A W(o) = lim e^lpf c (0"). 

€ — >OQ £ — 'OO 

Using the difference equation T(z + 1) = zT(z) it is straightforward to verify that 
lim c -e <+(* + 1 ) = I d^.+ jz + l) 

In addition, p BC = p B . It follows that 

Ag g ( P ^ + A) _ o|A | Ag^+A) 



lim e~ e|A| 2 2|A| t VP = 2 

A B C( p BC) Af(p B ) 



e — >oo 



which clearly implies the statement. □ 

We proceed to deduce the recurrence relations in question for the multivariable 
Bessel polynomials. As we will see below, these recurrence relations have a simple 
expression in terms of the functions 

~b, x K + z ~b i \ ((a - l)/2 + z) 

V (Z) = , W (Z) — ; r -. 

y ' z ' K ' 2z(2z+l) 

It is readily inferred from the difference equation T(z + 1) = zT(z) for the Gamma 
function that these functions appear as the coefficients in the following difference 
equations for the functions d B ± and d B ± : 

Lemma 6.1. The functions d B ' ± and d B ± satisfy the difference equations 

d B + (z + 1) = v B (z)d B + (z), d B _(z + 1) = v B (-z l)d B _(z), 

+ !) = d B t _(z + 1) = -w B (-z - l)d B t _(z). 

We are now ready to state and prove the simplest recurrence relation for the 
multivariable Bessel polynomials, corresponding to the elementary symmetric poly- 
nomial of degree one. 
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Proposition 6.2. The multivariable Bessel polynomials Y x satisfy the recurrence 
relation 

~(xi + ■■■ + x n )Y x (x) =^Vi (p B + A) (9 XW (x) - Y x (x)) 

(22) 

+ (^ + A ) (Yx it) (x)-Y x (x)) 

i 

with 

(23) V ±i (z) = w B (±z i )l[v B (±z i + z j )v B (±z i -z j ), 

where the first and second sum is over all i = 1, . . . , n such that \{ % > and \n\ is a 
partition, respectively. 

Proof. We fix a partition A = (Ai , . . . , A n ) . In order to make use of the orthogonality 
of the multivariable Bessel polynomials, as established in Theorem 17. 11 we assume 
that k > and that the condition (|2"T|) is satisfied for m = |A| + 1. It follows from 
(|3ej) and condition (1) in Definition 13.11 that 



^ n n 

(24) -( Xl + ■ ■ ■ + x n )Y x (x) = %Yxm {x) + %Y x (x) + £ V^% } (x) + P(x) 



= 1 i=l 



for some coefficients Vi and a symmetric polynomial P of degree less than |A| — 1. We 
let jj, = (/ii, . . . , fi n ) be a partition such that < |A| — 2. Then, (xi+- ■ ■ + x n )Y fl (x) 
is a linear combination of multivariable Bessel polynomials Y v (x) such that \v\ < |A|. 
It thus follows from Theorem 17 . 1 1 that 



<|A|-2 



= ((xi + ---+x n )Y x M =(p,y, 

Since the multivariable Bessel polynomials K, with \fj,\ < |A| — 2 span A 
(c.f. the proof of Corollary 15. ip . this implies that P = 0. Moreover, by setting 
x = (0™) in ([24]) and using the specialisation formula (f2Tj) wc find that 



71 

V = -^2(Vi + 



i=l 

Hence, the expansion of (x% + ■ ■ ■ + x n )Y x (x) in 'renormalised' Bessel polynomials 
Y" M is indeed of the form ([22]) . There remains to compute the coefficients V±i 
for i = 1, . . . ,n. It can be inferred from (6.24) in Chapter VI of Macdonald's 
book [Mac95] that 

{xx+'-- + x n )P x {x) = ip'xw / A P AW (x) 

i 

with 

,/ _ tt n(j - i + 1) + Aj - Xj n(j - i - 1) + A t - Xj + 1 
^aw/a - 11 K (j — i) + A* — A~ k(j - i) + A, - Xj + 1 ■ 

Applying this recurrence relation to the left hand side of l|22[). and then equating 
coefficients with the right hand side we obtain 



A^+aW) 
Vi ~ A B ( P * + X) ^> x 
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for i = 1, . . . ,n. Using Lemma l6.1[ and the formula for ^'xw /\ stated above, it is 

now straightforward to verify that the coefficients Vi indeed are given by (|2"5)) for 
i = 1, . . . , n. That (|23p is valid also for i = — 1, . . . , —n can be verified in a similar 
manner, using the series expansion of the multivariable Bessel polynomials obtained 
in Proposition ^. 21 This proves the statement under the assumed conditions on the 
parameters a and k. However, by analytic continuation in these parameters the 
proof immediately extends to any values of a and k for which the multivariable 
Bessel polynomial in question is well-defined. □ 



In principle we could use the same method as in the proof of Proposition 16 . 21 to 
deduce higher order recurrence relations for the multivariable Bessel polynomials. 
However, it is clear that this would lead to lengthy and complicated computa- 
tions. We will therefore use a different approach. To this end we recall that, by 
exploiting limit transitions from the multivariable Askey- Wilson polynomials, van 
Diejen [vD99] obtained recurrence relations for a number of different families of 
orthogonal polynomials. Among them are the BC n Jacobi polynomials. As we will 
show below, by applying the limit transition obtained in Proposition 14.11 to these 
latter recurrence relations it is simple and straightforward to obtain the remaining 
recurrence relations for the multivariable Bessel polynomials. We remark that van 
Diejen used a definition of the BC n Jacobi polynomials - in his notation denoted 

- which is somewhat different from the one employed in this paper. It is readily 
verified that the precise relation between these two definitions is given by 

V\{ixi, . . . , ix n ; v , v\, v) = Pjf c (x ll ...,x n ; h,k 2 , k 3 ) 

with the parameters [vq^vx^v) = (fci + k 2 , k 2 , k 3 ), and where i here denotes the 
imaginary unit. In order to avoid the use of more than one definition of the BC n 
Jacobi polynomials, we will formulate van Diejen's results in terms of the polyno- 
mials P\ ■ For that we require the 'renormalised' BC n Jacobi polynomials 

pBC _ q-2|A| Af C (p BC ) p BC 

It is clear from the discussion above that, as in the case of the multivariable Bessel 
polynomials, these 'renormalised' BC n Jacobi polynomials can be characterised by 
the fact that their constant term equals one. Furthermore, the functions 

_ ((fc 1 +2fc 2 )/2 + z)((fc 1 + l)/2 + z) 
V [Z) ~ z ' W (Z) ~ 2z(2z + l) 

and the symmetric trigonometric polynomials 

^ c (.) = (-ir+ i £ n sinh2 f> 

ic{i,...,n} iei 

\I\=r 

will be required. We note, in particular, that the trigonometric polynomials E r are, 
up to a difference in sign for odd r, the elementary symmetric polynomials in the 
functions ti(zi) = — sinh 2 # with i = 1, . . . ,n. Given an index set I C {1, ... , n} 
and a configuration of signs e» = ±1, i £ I, we will make use of the short-hand 
notation 
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For each r = 1, . . . , n van Diejen [vD99] (see Theorem 6.4) obtained the following 
reccurence relation for the BC n Jacobi polynomials: 

(25) E BC (z)P BC (z) = £ tW_ m (p BC + A) V eI IC {p BC + A) P BC ti Wj 

where the sum is over all index sets / C {1, . . . , n} with < |/| < r and configura- 
tion of signs Cj = ±1, i € /, such that A + e £ / is a partition, and where 



uj,m - {-i? e e (n 

KCJ c fc = ±l \keK 



w^(e k z k ) 

K 

\K\=p keK 

x II v BC (e k z k + e k ,z k ,)v BC (-e k z k - e k <z k > - 1) 

k,k'£K 
k<k' 

x j~[ v BC (e k z k + z i )w sc (e fe z fe - Zj) J , 

keK ) 
j£J\K 

Veij(z) = Y[w BC (eiZi) J[ v BC {e lZi + e v Zl ,)v BC ! (e 4 *j + e^/ + 1) 

i&I i,i'el 

x Y[ v 80 ^ + z^v^ieiZi - zj). 
je.J 

The recurrence relations for the multivariable Bessel polynomials will have the 
same structure as (|2"5"|) . We only have to make the appropriate substitutions. More 
precisely, with 

■l) p+1 



^=^- e n 

IC{l,—,n} i€l 
|J|=r 



we have the following: 



Theorem 6.1. For each r = 1, . . . , n the 'renormalised' multivariable Bessel poly- 
nomials Y\ satisfy a recurrence relation of the same form as (|25| . but with E^ 
substituted for E BC , p B for p BC , and the functions v B and w BC replaced by v B 
and w B , respectively. 

Proof. We let the values of the parameters (^1,^2,^3) be given by (fT4| . Then, 
i) BC = v B and p BC = p B . In addition, it is readily verified that 

lim e- f -w BC (z) = -w B (z), 

e—>oo 2 



lim e-^E BC ( Zl , ...,z n ) = ^E B {e z \. . . , e z "). 



1 

It follows from Proposition 14. II and the proof of Proposition 16.11 that 
lim P BC ( Zl +e,..., z n + e)=Y x {e z \...,e z "). 

e — >oo 

We thus obtain the statement by first multiplying both sides of ([23]) by 2 r e _er and 
then taking the limit e — > 00. □ 
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7. Orthogonality and norms 

In the one-variable case it is well known that there exist no (positive) measure 
on the real line with respect to which all Bessel polynomials are orthogonal; see 
e.g. Section 18 in Krall and Frink [KF49]. However, for a fixed maximum degree, 
and sufficiently negative values of a, the corresponding Bessel polynomials can be 
shown to be orthogonal with respect to just such an inner product. In this section 
we will establish a similar statement for the multivariablc Bessel polynomials. More 
precisely, under the assumption that k > 0, we will prove the following: given any 
non-negative integer m, and sufficiently negative a, the Bessel polynomials Y\ of 
degree at most m form an orthogonal system with respect to the inner product in 
the Hilbert space L 2 (W%, d^) with 

n 

d^ B \x; a, k) = Y[ ir 2 e" 2/l - n \xi - x 3 \ 2K dxi ■■■dx n - 

i—l i<j 

Suppose now that we apply the Gram-Schmidt orthogonalisation procedure to this 
inner product to obtain a set of symmetric polynomials labelled by the partitions 
A = (Ai, . . . , A n ) of weight |A| < to. It is then clear that, with the appropriate 
normalisation, these symmetric polynomials will coincide with the multivariable 
Bessel polynomials for the parameter values mentioned above. Moreover, since the 
multivariable Bessel polynomials are rational functions in a and k (see Corollary 
13. 1|) this fact can be immediately extended to generic parameter values. 

At this point we fix a non-negative integer m, and consider the resulting subspace 
A^™ 1 C A„, as defined in Section 2.1. We let (•, -) atK denote the inner product in 
£ 2 (R^,d^ B )), i.e., 

(26) (f, g ) atK = [ f{x)g\x)d^ B \x-a,K) 

for any f,g G L 2 (R" , d/i^). It is clear that whether or not A^ m is contained 
in this Hilbert space depends on the values of the parameters a and k. A precise 
condition is given by the following: 

Lemma 7.1. Assume that k > 0. Then A| m C L 2 (K™ , d^) if and only if 

(27) a < -2(m + re(n - 1)) + 1. 

Proof. We fix two partitions A = (Ai, . . . , A„) and /i = [fix, . . . , /x„) such that 
|A|, \(t\ < rn, and consider the integral 

(28) f mxm„d^ B K 

Since d/j>( B \a, k; x) is symmetric in x, we can replace the domain of integration by 
the open subset in M™ defined by the inequality 

< x n < x„_i < • • • < x\. 

We observe that 

2k 

dx\ • • ■ dx n . 



d^ B \a,K ] x) = \{\ 



-2+2 K (n-j) e -^- 



n 



1_ Zi 

J J j 
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0<a:„<"'<xi 



It follows that 

mx{x)m^{x)d^S B ' 1 (a, k; x) 

r - 2 

< / m x (x)m, l {x)Y[\x i \ a ' 2+2K( - n ^e--dx 1 ---dx n . 

Ji)<i„<-<ii i=1 

Clearly, the latter integral exists if and only if 

Ai + /ii + a - 2 + 2n(n - 1) < -1. 

By setting A = fj, = (m) we thus conclude that A^ m C L 2 (M™ , d/i^) if a satisfies 
([27]) . We now suppose that a does not satisfy the condition (f27|) . Wc keep A = \i = 
(m), and fix all variables x except X\. It is then readily seen that the integrand of 
can not be integrated in x\. Consequently, the integral does not exist. □ 



In order to establish that the multivariable Bessel polynomials contained in A^ m 
form an orthogonal system we proceed to prove that their eigenoperators, obtained 
in Section 5, are symmetric operators on A^ m with respect to the inner product 

("l ')a,K- 

Lemma 7.2. Assume that k>0 and that a satisfies the condition (|27p . Then 

( D *f>9) a , R = {f,D%g) aR 

for all f,g € A^ m and d = 1, ... ,n. 
Proof. We let 

n 

W(x;a, K ) = l[xr 2 e- 2 ^l[\ 



Xl X 1 



°3\ 

i=l i<j 

so that d^ B \a, n;x) = W(x; a, n)dx\ ■ ■ ■ dx n . We will prove the statement by 
verifying that 

n n 

(29) E E <^V> Kei9) a , K = E E <^i.-/. D ^i9) a , R 

e=±l i=l e=±l i=l 

for all /,j e AJ==™ 1 , and positive and non- negative integers p and q, respectively. By 
partial integration in Xi we deduce from (|26p that 

x iJT~Dp_ li f, Dq^g 



a , k 



9a:.. 



It is readily verified that 



W - x -jL Xi W = a - 1 + — + «(n - 1) +K ^^, 

OXi Xi . . Xj 
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The definition, as stated in (fl6|) . of the differential operators thus imply that 
{D* ei f,Dl ti g) aK = - (DB_ l €i f,ex~Dl ei g 

(D*_ liei + Dp_ lt _ ei )f, |(a-l + A) D^g 

^((D^+D^f^^D^g 



2 



a . k 



Taking the sum over e and i, and using the invariance under the interchange of the 
summation indices i and j, as well as under the substitution of — e for e, we obtain 
(123. □ 



To recapitulate, we know from Theorem l5.1l that the multivariable Bessel polyno- 
mials are common eigenfunctions of the differential operators for d = 1, . . . , n. 
In addition, we have just shown that these operators, for appropriate values of k 
and a, are symmetric on A^ m with respect to the inner product (•, -) a , K - Since the 
differential operators in the algebra 2$ n = C[Z?f , . . . Df % ] separate the multivari- 
able Bessel polynomials in question, as shown in Corollarv l5.31 we thus obtain the 
following orthogonality result: 



Theorem 7.1. Assume that n > and that a satisfies the condition (J27J) . Then the 
multivariable Bessel polynomials Y\ with |A| < m form an orthogonal system with 
respect to the inner product (•, -) a ,K, for any two partitions A = (Ai, . . . , A„) 
and \i = . . . , /i„) such that |A|, \n\ < m, 

(Y x ,Y li ) a , lt = 0, ifX^fi. 

We proceed to compute the squared norms 

N x (a,n) = (Y x ,Yx) a , K , 

thus enabling us to convert the multivariable Bessel polynomials contained in A^™ 
into an orthonormal system. We will proceed in two steps: first we use the recur- 
rence relations obtained in Section 6 to reduce the problem to that of computing 
A^q), this latter normalisation factor is then computed using a well-known limit- 
ing case of an integral formula due Selberg [Sel44]; see also Forrester and War- 
naar [FW07]. The first step is accomplished in the following: 

Lemma 7.3. Assume that n > and that a satisfies the condition (|27p . Let 
A = (Ai, . . . , A„) be a partition such that |A| < m. Then 

m * ( lV A| Af(^ + A)A^ + A) 

Proof. We will prove the following equivalent statement: the squared norms 
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of the 'renormalised' multivariable Bessel polynomials Y^ with < m, as defined 
in (f20|) . are given by 



We note that each partition of length at most n can be written as a unique linear 
combination with positive integer coefficients of the integer vectors u> r = ei+- • -+e r 
with r = 1, . . . , n. To prove the statement it is thus sufficient to verify that both 
sides of (|31[) transform in the same manner under the replacement of p, by p + c<j r 
for each r = 1, . . . ,n. Using Theorems 16. II and 17. II we deduce that 



N M+Ur 



^{t...,r},{r + l,..,r l }(/ ?B +/ i ) 
{l,...,r},{r+l,...,n} 



^{l,..,r},{r+l,...,«}(-P B 'I 1 ' U J f, 



^ ■ i'i 

where we used the expansion of E B Y„ to obtain the first equality, and that of 
E B Y^ +UJr to obtain the second. It follows from Lemma HOI that 

yB (z) Kjz+UJr) 

V {l,...,r},{r+^...M^ Z > ~ Jfi^ ' 

and that 

A B (z + uj r ) 

We thus conclude that 

M A B (p B + » + LJ r ) = V { B l,...,r},{r+l,...,n}(-p B -V-^r) A B (p B +p) 

A B (p B +p + co r ) V" { 1..,, } , {J . +1 ,...,„ } (P B +M) k B (p B + fi)' 

and the statement follows. □ 

There remains only to compute the normalisation factor . 
Lemma 7.4. Assume that k > and that a < —2n(n — 1) + 1. Then 

T(iK)T(-a - n(n + i - 2) + 1) 



(32) N {0) = 2 «"("- 1 )+( a - 1 )"n! ]J • 



r(«) 



Proof. We first observe that Lemma I7TT1 ensures that 7V( ) is well-defined. Changing 
variables to (y\, . . . , y n ) = (2/x\, . . . , 2/x„) we find that 

poo poo n 

jy (0)=2 '"M)+M)»/ d2/1 .../ % l n^ 2K( ^ i) ^ni^-%i 2K - 

- 70 J ° i=l i<j 

As indicated above, this integral is a well-known limiting case of the Selberg integral. 
Indeed, by specialising to the root system B n in Conjecture 6.1 in Macdonald 
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[Mac82] it is readily inferred that 

/>oo />oo n 

/ '///,••• / </.</, IK ' " 1 < "II th 1 " 

Jo Jo -_i 

n 

from which the statement clearly follows. □ 

We proceed to combine Lemmas 17.31 and \7A\ to deduce a fully explicit expression 
for the squared norms N\. 



T(in)T(a- 1)) 



Proposition 7.1. Assume that k > and f/iai a satisfies the condition (127|) . Let 
A = (Ai, . . . , A„) 6e a partition such that |A| < m. Then 



N x = c x 11 

l<i<j<n 



X 



r(«(j - t + 1) + Aj - Aj)r(«(j — i — 1) + 1 + Aj — Aj) 
r(/s(j - *) + Ai - AjOrf/sO' - *) + 1 + ^ - Aj) 

T(— a — ft(2n — i — j) + 2— Xi — Xj) 
X r(-a - rc(2ra - i - j + 1) + 2 - A, - A^) 

r(— a — n(2n — i — j) + 1 — A 2 - — Xj) 
X r(-a - K(2n — i — j — 1) + 1 — Aj — Aj) 

^- / r(-a - 2k(« - «) + 2 - 2A t )r(-a - 2k(ti - i) + 1 - 2A») 



n 



. x r(-o - «(n - i) + 2 - A.,) 



x r(«(n -i) + 1 + Ai) 



with 



Ca — 2 2 l A l+ KrI ( n_1 )+( a_1 ) n 77,! 

In particular, it is manifest that N\(a, k) is a positive analytic function of a and k 
(under the assumed conditions on these parameters). 

Proof. We will use the notation (a) m for the falling factorial, defined by (a)o = 1 
and 

T(a + 1) 

(a) m = a{a - 1) ■ ■ ■ (a - m + 1) = — — — 

1 (a — m + 1) 

for m > 0. In addition, we let (ai, . . . , a p ) m — (ai) TO • • • (a p ) m , and similarly for 
the Pochhammer symbol [a] m . It follows from the formulae T(a + m)/T(a) = [a] m 
and (— l) m [a] m = (— a) m that 



( i)M Af(p B + X)A B ( P B + X) = j-j [<J * + 1), «(* - i - 1) + l] A< _ Aj 



A? (p b )A s (^) ^ [K(j i), n{ 3 - i) + 1] As _ Aj 

( — a — fc(2n — i — j + 1) + 1, —a — «;(2n — i — j — 1)) A . +A . 
( — a — k(2h — i — j) + 1, —a — n(2n — i — j)) A+A . 

[n(n — i) + l] A ( — a — «(n — i) + l) ; 



n 



( — a — 2k(u — i) + 1, —a — 2ft(n — «)) 



2A, 
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Using the identity 

n r(a(j-i + l) + /3) _ -A- T(±ia + 0) 
iM<n r H3-i) + 0) l\ r((l-i±i)a + /J)' 

valid for generic a and it is readily verified that 

tt r(«(j - z))r(«(j -») + i) A r(«)r(«(i - 1) + 1) 

l<i< 3 <n r(K{ 3 - 1 + - 1 - ^ + V tl ^ 

and that 

yr r(-a - k(2u — i — j + + 2)T(-a - n(2n -i-j-V) + l) _ 
Ki^ <n r(-a - k{2h -i-j) + 2)T{-a - n{2n - i - j) + 1) ~ 

-q r(-a - K(2n - i) + 2)r(-o - K(2n - 2i) + f ) 



r(-a - n{2n - 2i + 1) + 2)r(-a - «(2n - i - 1) + 1) ' 

By cancelling common factors in the resulting nominator and denominator, and 
using the difference equation T(z + 1) = zT(z), it is now a straightforward exercise 
to verify that Lemmas 17.31 and 17.41 combine to give the statement. □ 

8. Concluding remarks 

In this final section we briefly return to the discussion of the relation between 
the multivariable Bessel polynomials and the hyperbolic Sutherland model with 
external Morse potential, and present some remarks on the problem of constructing 
a (moment) functional with respect to which all multivariable Bessel polynomials 
would be orthogonal. 

8.1. Eigenfunctions of the hyperbolic Sutherland model. We will through- 
out this section assume that k > 3/2. To each permutation P G S n we associate 
the following open subset of K™: 

A P = {z e W 1 : zp ( i) < ■ • • < Zp(„)}. 

In addition, we assume that the parameters a and k are such that there exist at 
least one partition A = (Ai, . . . , A n ) which satisfies the inequality 

(33) o<-2(|A| + K(n-l)) + l. 

For each such partition A and permutation P G S n we have the following eigen- 
function of the Schrodinger type operator |(2| : 

(34) ^\ z) = {fM^ Y ^ *eA P> 



0, z^Ap, 

with N\ as in Proposition 17. II It is clear from the discussion in Section 7 that this 
eigenfunction is normalised to one. We remark that all the second order derivatives 
of the eigenfunction (|3"4"|) are contained in the Hilbert space L 2 (R ra , dzi ■ ■ ■ dz n ) 
precisely for k > 3/2. In addition, it seems that for k < 3/2, or rather for 2k(k — 
1) < 3/2, it is a rather delicate matter to associate a domain to such that 
the resulting Schrodinger operator is self-adjoint; c.f. Feher et al. [FTF05] which 
contains a detailed investigation of the equivalent problem in the case of the three 
particle Calogero model. 
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The eigenf unction (j34|) of the Schrodinger type operator ([2]) has the following 
natural physical interpretation: it represents a bound state in the corresponding 
quantum many-body system where all particles are confined to the region Ap. 
Indeed, the fact that 'J'o vanish at each hyperplane Zi = Zj with i, j — 1, . . . , n such 
that i j implies that the probability flux across such a hyperplane is zero; see 
e.g. Section 2.4 in Sakurai [Sak94] for a definition of the probability flux. In more 
concrete terms this means that the particles can not pass each other. Consequently, 
once confined in a given region A p they will remain there indefinitely. For a further 
discussion of this point see e.g. Calogero [Cal71]. 

We note that for given values of the parameters a and k, there arc only a finite 
number of partitions A that satisfy (|33|) . i.e., the Schrodinger type operator ((2]) has 
only a finite number of eigenfunctions of the form (|34|) . This is a manifestation of 
the fact that the Morse potential only supports finitely many bound states. 

8.2. The multivariable Bessel moment problem. As we mentioned at the 
beginning of Section 7, in the one-variable case it is well known that there exist no 
(positive) measure on the real line with respect to which all Bessel polynomials are 
orthogonal. This naturally leads to the question of whether they form an orthogonal 
system in some more general sense. To be more precise, one might inquire whether 
or not there exists a non-trivial (moment) functional 

/ : A„ -> C 

with respect to which all multivariable Bessel polynomials are orthogonal. For 
one variable this problem is known to have a positive solution, and a number of 
explicit integral representations of such a functional have been constructed. We 
mention, in particular, the work of Krall and Frink [KF49] , who obtained an integral 
representation with integration along the unit circle, and Evans et al. [EEKKL92], 
in which the integration is along the positive real line. In this section we will briefly 
consider this problem in the two-variable case. 

We thus set n — 2, and recall that the algebra A2 is freely generated by the 
elementary symmetric functions 

ex = xi + x 2 , e 2 = x\x 2l 

see e.g. Section 1.2 in Macdonald [Mac95]. It is readily verified that, as a linear 
operator on A 2 , D B is given by 

D B = (e\ - 2e 2 )dl + 2 ei e 2 d ei d e2 + 2e 2 2 d 2 e2 

+ (1 + (2k + a)ei)d ei +{e 1 + 2{n + a)e 2 )d e2 , 

where we have used the notation d ei = d/dei, and similarly for d e2 . For a functional 
I : A2 — > C and p € A 2 we define d ei I, d e2 1 and pi by setting 

9 ei I(q) = -I(d ei q), d e2 I(q) = -I(d e2 q), pl(q) = I(pq), 

for all q 6 A 2 . With this definition in place, it is a matter of a straightforward 
computation to verify that DB is symmetric with respect to /, i.e., I(qD B p) — 
I(pD B q) for all p, q S A 2 , if and only if / satisfies the differential equations 

d ei (el - 2e 2 )I + d e2 e x e 2 I - (1 + (2k + a)e x )l = 0, 
d ei e x e 2 I + 2d e2 e\l - (e x + 2(k + a)e 2 )I = 0. 



MULTIVARIABLE BESSEL POLYNOMIALS 



27 



We let I nrn = / (e^e™) be the moments of I with respect to the linear basis of A2 
formed by the symmetric polynomials e™e™ with n,m £ N. The differential equa- 
tions for / stated above then translate into the following two recurrence relations 
for the moments I nm ' 

(35a) (n + m - (2k + a))I n+ i tTn - 2nl n ^ liin+ i - I nm = 0, 

(35b) (n + 2m - 2(k + a))I n . m+1 - J n+ i, m = 0. 

By induction in n + m it is a simple exercise to show that once the moment /o,o is 
fixed these recurrence relations uniquely determine the remaining moments I nm . It 
is easily seen that for generic values of the parameters a and n (and n = 2) the eigen- 
values of D B are non-degenerate. Assuming this to be the case, the two-variable 
Bessel polynomials thus form an orthogonal system with respect to a functional I 
determined by the recurrence relations (I35|) and the moment To.o- 

From the discussion above it is not clear whether such a functional / might 
have a natural integral representation, or even if its moments (with respect to the 
elementary symmetric functions) have a simple expression. It is furthermore not at 
all clear what happens for n greater than two. We plan to return to these questions 
in a forthcoming paper [Hal] . 
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